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The Hammersley process relates to the statistical properties of the maximum length 
of all up/right paths connecting random points of a given density in the unit square 
i from (0,0) to (1,1). This process can also be interpreted in terms of the height of 

! the polynuclear growth model, or the length of the longest increasing subsequence 

, in a random permutation. The cumulative distribution of the longest path length 

Q , can be written in terms of an average over the unitary group. Versions of the Ham- 

. mersley process in which the points are constrained to have certain symmetries of 

. the square allow similar formulas. The derivation of these formulas is reviewed. 

Generalizing the original model to have point sources along two boundaries of the 
^ . square, and appropriately scaling the parameters gives a model in the KPZ univer- 

. sality class. Following works of Baik and Rains, and Prahofer and Spohn, we review 

. the calculation of the scaled cumulative distribution, in which a particular Painleve 

■ II transcendent plays a prominent role. 

m : 
o ■ 

^ ! 1 Introduction 

■ The aim of this review is to explain aspects of developments over the past few years relating 
' some observables in statistical mechanics models to random matrix averages and then to Painleve 

S ■ transcendents. In addition to the theoretical interest in these inter-relationships, the fact that the 

' Painleve transcendents are readily computable means that quantitative predictions for certain 

^ ' order one scaled observables are available for the first time. Perhaps the most significant such 

d ' result is the calculation due to Prahofer and Spohn [33] of the exact two-point scaling function 

for one-dimensional stationary KPZ (Kardar-Parisi-Zhang) growth. KPZ growth is generally 
believed (see e.g. |251 |Hj) to underlie a diverse number of growth models in 1 + 1 dimension. 
One model within the KPZ universality is the polynuclear growth (PNG) model. The particular 
variant of the latter relevant in this context was solved in terms of a Painleve transcendent 
known from random matrix theory |27j by Baik and Rains [7j, and it is this solution which is 
interpreted and computed in [521 • Section 4 of the present work gives some details of the exact 
solution. 

We begin in Section 2 by reviewing the calculation of the cumulative distribution for the 
longest path in the Hammersley process. This will be shown to be equivalent to computing 
the cumulative distribution for the maximum height in the PNG model, the longest increasing 
subsequence length for a random permutation, or the maximum displacement of certain families 
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of non- intersecting paths. The cumulative distribution is given as a particular random matrix 
average over the unitary group. In Section 3, four different symmetrizations of the Hammersley 
process are considered. The cumulative distributions of the longest path in each of these cases can 
again be written as particular random matrix averages, involving the orthogonal and symplectic 
groups in two of the cases, and the unitary group in the remaining two. After presenting some 
details of the calculation of the scaled distribution for the variant of the PNG model of relevance 
to KPZ growth in Section 4, we conclude in Section 5 by indicating aspects of the Painleve 
transcendent content of the averages over the orthogonal and symplectic groups encountered in 
Section 2. 

2 The Hammersley process 
2.1 Relationship to permutations 

The Hammersley process (see E] for an extended account of different emphasis to that given 
here, and for references to the original literature) refers to the following stochastic model. In the 
unit square mark in points uniformly at random according to a Poisson rate with intensity A^, 
so that the probability the square contains points is equal to X"^^ e~'^^ /N\. Form a continuous 
path by joining points with straight line segments of positive slope, which are thus orientated 
up and to the right. Extend this path to begin at (0, 0) and finish at (1, 1) by adding an extra 
segment at both ends, and define the length of the extended path as the number of points it 
contains. Take as the primary observable quantity the stochastic variable, say, specifying the 
maximum of the lengths of all possible extended paths (see Figure ^ . 
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Figure 1: Eight points in the unit square, and the extended directed paths of maximum length. 
Since the number of points in these paths equals three, here = 3. 

For any particular realization of exactly N points the Hammersley process gives a geometrical 
construction of a random permutation of {1, 2, . . . , N}. This comes about by first labelling the 
X coordinates of the points byO<xi<---<j;Ar<l and similarly the y coordinates by 
< yi < • • • < Vn < 1- Each point will then have a coordinate of the form {xj,yp(^j-j) where 
{P(l), . . . , P{N)} is a permutation of {1, 2, . . . , N}. The quantity also has an interpretation 
in terms of the permutation. Thus the analogue of an up/right path connecting points is a 
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subsequence 1 < ii < i2 < • • • < jV < such that -P(ji) < -P(i2) < ••• < P{jr)-, which is 
referred to as an increasing subsequence. The length of an increasing subsequence is defined 
as the value r. We then see from the definitions that the maximum length of all increasing 
subsequences of P coincides with 

2.2 Polynuclear growth model 

Consider the x-t half plane t > 0. Let this half plane be filled with points uniformly at random 
and such that the mean density is unity. These points are to be thought of as seeds for nucleation 
events of layered growth. In a droplet model, at (x, t) = (0, 0) a single layer, taken to have zero 
height, starts spreading with unit velocity to the left and to the right. Forming on top of the 
ground layer are new layers of unit height. These layers, or parts thereof, are formed at space- 
time positions {xi,ti) for each nucleation event bounded by the 'lightcone' axis u = {t + x)/^/2, 
V = (t — x)/ \/2; nucleation events outside this cone are not created at a time that their position 
coordinate makes contact with the ground layer or its growth. The nucleation events {xi,ti) 
inside the lightcone create the beginning of a portion of a layer of unit height on top of the 
ground layer, or existing layers, at position Xj. The layers are formed by the growth of the 
nucleation events with unit velocity to the left and to the right; if two growing portions of a 
layer collide, then growth at that point ceases and the two portions become one, growing only 
at the end points of this one portion (see Figure [21 for an example). Of interest is the statistical 
properties of the height at the origin after this growth process — known as the polynuclear 
growth (PNG) model — has been underway for time t = T. 



Figure 2: Example of the plateau profile at the time of four successive nucleation events, 
including the initial event (which is labelled the 0th event and its plateau the 0th level). Note 
that between the second and third nucleation event, two plateaus on the first level have coalesced. 

The first observation is that only those nucleation events in the region [u = 0,u = T/V2] X 
[v = 0,v = T/\/2] of the lightcone can contribute to the height at x = up to time t = T. 
Suppose in a realization of the nucleation events there are N points in this region. For a Poisson 
process of unit density this occurs with probability e~^^ /N\, where = is the area of 
the region. Use the construction of the previous subsection to associate with the configuration 
of points a permutation P (see Figure ^ . Also indicated in Figure Ol are the world lines of 
the nucleation events, which show clearly the layered structure of the growth, and in particular 
the height at the origin after time T. The layers in which the particular nucleation events 
occur are simply related to the permutation P. This is done by partitioning the permutation 
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into decreasing subsequences using the leftmost digits at all times. The jth such decreasing 
subsequence corresponds to the jth layer in the growth process. For example, in Figure 01 the 
permutation is 5374162, and the decreasing subsequences formed from the leftmost digits are 
(531)(742)(6). It is easy to see that in general the number of decreasing subsequences of this 
type is equal to the length of the longest increasing subsequence of the same permutation. Thus 
studying the height at the origin in the PNG model after time T is equivalent to studying the 
maximum path length in the Hammersley process with intensity = 



t 




Figure 3: World lines for the endpoints of the plateaux. Only nucleation events inside the 
square shaped region including the lines x = it, and the lines from t = T to these lines affect 
the height at the origin. The nucleation points occur at v shaped configurations, while the 
inverted v part of the worldlines correspond to the joining of the plateaux originating from two 
different nucleation events. The labelling on the lines x = t and x = —t allow the world lines to 
be identified uniquely with a permutation. 

2.3 Robinson-Schensted-Knuth correspondence 

To analyze the Hammersley process requires a formula for the number of up/right paths less 
than some prescribed value / say. For this purpose, in a realization containing A'^ points, we first 
associate with each point (xj,yp(j)), j = 1,. . . ,N the permutation matrix defined so that the 
entries {j,P{j)) are equal to unity with all other entries equal to zero. We then apply what is 
essentially Viennot's shadow method 40^ to give a bijection between permutation matrices and 
certain configurations of lattice paths, the outermost member of which can be interpreted as the 
profile of a lattice variant of the PNG model jSO] . 

In fact it is possible to give a bijection between a n x n non-negative integer matrix X = 
[xij]ij=i,..,,n and a pair of so called up/ right horizontal (u/rh) non-intersecting lattice paths. 
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This bijection is equivalent to the celebrated Robinson-Schensted-Knuth correspondence map- 
ping X to a pair of semi-standard tableaux The non-intersecting lattice paths are defined on 
the square lattice and start at x = 0, one unit apart in the y-direction at y = 0, . . . , — (n — 1), and 
finish at = n — 1, with y-coordinates fii — {I — 1) (/ = 1, . . . , n) where /^i > /i2 > • • • > fJ-N ^ 0. 
The path starting at y = — (/ — 1) is referred to as the level-/ path. Each path may move either 
up or to the right along the edges of the lattice, with the constraint that the paths may not 
intersect. When X is a permutation matrix and so has exactly one non-zero entry, equal to 
unity, in each row and column the paths are further restricted so that for each allowed x-value 
there is exactly one up step in total, which is of unit length. 

In the general case the entries Xij of X, where for convenience the rows are labelled from 
the bottom, represent the heights of columns of unit length centred about x = j — i which occur 
at time t = i + j — This labelling is simple to implement by first rotating the matrix 45° anti- 
clockwise. The columns are to be placed on top of the level-1 path formed by earlier nucleation 
events and their growth. On this latter point, during each time interval the existing profile 
or profiles are required to grow one unit to the left or to the right, with any resulting overlap 
recorded on the path one level below. What results are a set of at most n non-intersecting 
lattice paths, equivalent to a pair of u/rh non- intersecting lattice paths each with the same 
final positions. Furthermore, it is easy to see that the process is invertible, and so there is a 
bijection between non-negative integer matrices and pairs of u/rh non-intersecting lattice paths. 
As already remarked, in the case of permutation matrices, these paths have the additional 
constraint of containing exactly one up step in total for each allowed x-value, x = ±{2n+3/2—2j) 
{j = 1 , . . . , n) . An example is given in Figure |^ 
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Figure 4: The RSK mapping in the non-intersecting paths picture from a permutation matrix 
to a pair of u/rh non-intersecting paths, constrained so that there is exactly one up set in total 
for each allowed value of x. 

Non- intersecting u/rh lattice paths, constrained to have exactly one up step in total for each 
allowed x-value can be encoded as standard tableaux. The latter consist of an array of unit 
boxes stacked across rows of length //i, /X2, . . . , /^n respectively with the first box in each row 
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contained in the first column of the array etc. The standard tableau is said to have shape fj,. 
One plus the value of the x-coordinates that contain up steps are marked in order along row /. 
Thus the numbers in the array have the property of being strictly increasing both along rows 
and down columns, with each number from {1,2, . . . ,n} recorded once. The standard tableau 
is said to have content n. 

Let denote the number of standard tableaux of shape fi and content N. Let denote 
the longest increasing subsequence length in a realization of the Hammersley process containing 
N points. Then it follows from the correspondence with pairs of constrained u/rh lattice paths 
and thus pairs of standard tableau that 



(2.1) 



As a first step towards evaluating the sum in ()2.1() . note that the symmetry between rows and 
columns in a standard tableau implies we can write 



(2.2) 



:u'<i 



where /x'^ denotes the length of the first column. In terms of lattice paths, this can be understood 
from the bijection between u/rh lattice paths and dual u/rh lattice paths, where if in the former 
the maximum displacement is I, then in the latter there are exactly I lattice paths (see Figure 
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Figure 5: Drawn in heavy lines on the square lattice is a family of two u/rh lattice paths, 
constrained so that there is exactly one up step in total for each x-value, and with maximum 
displacement 3 units, while drawn in heavy dashed lines are the corresponding three dual lattice 
paths. Also given is the standard tableau encoding the two original u/rh lattice paths. 



2.4 Relationship to the lock step model of vicious walkers 

The task of computing the number of pairs of constrained u/rh lattice paths, each with the 
same final positions, and containing exactly / lattice paths can be considered as a special case 
of a more general counting problem. First we remark that a set of /* < / lattice paths between 
x = 0to3; = n — 1 can be uniquely converted to a set of / lattice paths by drawing horizontal 
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segments from x = Otoa; = n — 1 along y = —(/* — 1), ...,—(/ — 1), so to compute 1)2. 2|) the 
case of precisely I lattice paths is what is relevant. We generalize the rules of constructing the 
lattice paths so that for each allowed x value a; = 0,l,...,n — 1 there is again only one segment 
which is not horizontal, but now this segment may be either up or down, subject again to the 
non-intersecting constraint. These paths can be considered as the world lines for the stochastic 
evolution of / random walkers on a one-dimensional lattice. At each tick of the clock exactly one 
walker may move to the left or to the right one lattice space, subject to the constraint that this 
lattice site is not already occupied. This defines what is known as the random turns model of 
vicious walkers ^1 El ■ Let an up (down) step at a; = j — 1 be weighted Wj' {w~ ) , and suppose 
furthermore that each level-A; path returns to y = —{k — 1) at x = n — 1. Then with 

I 

Gi^n '■= (weight of level-/c path) 

paths k=l 

it is shown in that 

Gi,n = r de,... r n + v''^) n \^''' - ^''\"- (2-3) 

\ZT^) I. .^^ l<j<k<l 

But 

(/^)^ = G';^2Af wp=o,w+=i(p=i....,N) 

^i:^l\<l u)+=0, JDp =1 {p=iV+l,...,2JV) 

so we see from (|2.2p and (|2.3p that 

■ ^ ' ' ■'^'^ ■'^'^ j=\ l<j<k<l 

where {■)u(^i) denotes an average over the eigenvalue probability density function for random 
matrices from the classical group U{1), 

^<j<k<l 

To see the validity of the second equality in (|2.4j) . one first notes 

j=l p=0 ^ j=l j=l 

and then observes that only the p = N term in this sum is non-zero after averaging over U{1). 
Now, according to the definitions 



\27V 

Pr(F <l) = e-'' -W^'(^^ ^ ^) 

N=0 



m 



so substituting ([2^ gives [23 EH 



Pr(F</)=e-"Ye^Tr(j/+t/t)\ ^2.5) 
\ /ueu{i) 
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2.5 Relationship to Schur polynomials 

Consider the set of all semi-standard tableaux of shape A and content N, with each occurence 
of j in the numbering therein weighted wj. The total weight 

/ N ■\r^ #l's #2's #N's /r, r\ 

sx[wi, . . . ,wn) ■■= 2^ wf ■■■Wn (2-6) 

semi- standard tableux 
shape A, content N 

is a symmetric polynomial in wi, ... ,11)^ known as the Schur polynomial. Because in a standard 
tableaux each number occurs exactly once, it follows immediately that 

[wiW2 ■ ■ ■ wn]sx{wi, Wn) = Jn (2-7) 

where [wiW2 • • • wjy] denotes the coefficient of wiW2 • • • wn- This fact allows expressions such as 
involving to be evaluated as special cases of Schur function identities. Explicitly, in 
relation to (|2.1[) use can be made of the Schur function identity 

Af I 

s^iau aN)s^{h, . . . , 6iv) = ( n 11(1 + a,e*^'=)(l + 6, 6"*^'=))^^^^. (2.8) 

li:fj,l<l j=lk=l 

As to be revised in Section 4, this has direct relevance to the Johansson model. Our present 
interest is that by extracting the coefficient of ai • • • ajybi • • • 6jv from both sides, making use of 
()2.7|) on the left hand side, it follows that 

fj,:^i<l k=l 

in accordance with the evaluation of (|2.1j) implied by the first equality in (|2.4|) . 

2.6 A relationship to eigenvalue distributions 

It is of interest to note that the random matrix average over U{1) in (|2.5() also arises in an- 
other probabilistic setting jl4j . Thus consider the Laguerre unitary ensemble, specified by the 
eigenvalue probability density function proportional to 

N 

JjAfe-^' n Xi>0. 

1=1 l<j<k<N 

For a = ni — n2, ni > n2, it is realized by the eigenvalues of the random matrix X'^X, where X 
is an ni x n2 complex Gaussian matrix. Let £"2 (s;a;iV) denote the probability that there are 
no eigenvalues in the interval (0, s). Then it was shown in that the scaled probability 

A 



Er'iX;a):=lini^E^{-;a;N) (2. 
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is given by the right hand side of H2.5() with / = a and A ^ \/A (the superscript "hard" is used 
to denote the fact the eigenvalue density is strictly zero to the left of A = 0). This fact has been 
used in to give a straightforward proof that 

i^PK^^^^)=^^UE(.). (2.10) 

where -Fgue('S) is the scaled cumulative distribution of the largest eigenvalue for large random 
Hermitian matrices with complex Gaussian matrices, a celebrated result due originally to Baik, 
Deift and Johansson |Sj , which gives as a corollary the scaled distribution of the longest increasing 
subsequence of a random permutation. 



3 Symmetrizations of the Hammersley process 

3.1 Four symmetries of the square 

Baik and Rains [Zj have formulated and analyzed four symmetrized versions of the Hammersley 
process, in which the points are constrained to have particular reflection symmetries of the 
square. In these symmetrizations, the points are constrained to be invariant under reflections 
about the diagonal (0, 0) to (1, 1); diagonal (0, 1) to (1, 0); both these diagonals; and about the 
centre point (1/2, 1/2). The four cases are denoted 0, s, □ and h respectively. The first three of 
these cases can be further generalized to allow independent Poisson rates for points forming on 
the diagonal (s). Let us consider each case separately. 

3.2 The symmetry 

First we specify the Poisson process by which the points are added to the diagonal and below 
the diagonal. We start with a time interval [0, z], which is broken up into M smaller intervals 
of equal size, the latter being labelled j = 1, . . . , M. In each of these smaller intervals, add 
one point below the diagonal (together with its image above the diagonal) with probability 
z^(j — l/2)/M^ and a point on the diagonal with probability az/M. The probability that there 
are n points in the square is then given by the coefficient of in 

a aw\(^ (i- 1/2)^2 (j - l/2)z2u;2 



TT 1-- + — 1-^'"^/"^" 
iJ- V M )\ M2 ^ M2 

i=i 

It follows by taking the limit M ^ 00 in this expression that the probability of their being 
exactly n points after time z, when the points below the diagonal are added with rate zdz and 
those below the diagonal with rate adz is given by 

,-».,-.V2£:'f (3.1) 

n\ ^-^ \lmJ Z^ni\ 

m=0 

Furthermore, one sees that this probability conditioned so that there are n — 2m points on the 
diagonal is equal to the mth term in this sum. 
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Let denote the longest path length in a realization of the symmetrized Hammersley 
process consisting of a total of n points, m of which are below the diagonal. It follows from (|3.1j) 
and the sentence below that if the points are chosen according to the Poisson process specified 
above, then after time z 

00 „ [n/2] 

Pr(/^ <l)= e— E ^ E «"-'"^^n,™Pr(/^,„ < I). (3.2) 

n=0 m=0 

To compute Pr(/^^ < I) we first note that the sought realizations of the symmetrized Ham- 
mersley process are in correspondence with n x n permutation matrices constrained to have 
n — 2m non-zero entries on the diagonal, and which are furthermore symmetric about the diag- 
onal (recall that our convention is to count rows from the bottom and so the diagonal runs from 
the bottom left to the top right). It is straightforward to show that there are Sn,m distinct such 
permutation matrices. Next we make a correspondence between such permutation matrices and 
pairs of suitably constrained u/rh lattice paths, or equivalently pairs of suitably constrained 
standard tableaux. 

Now it is immediate from the rules of the PNG model that if the non-negative integer 
matrix X = [xij]ij=i^,„^n maps to a pair of u/rh paths {Pi,P2), then the transposed matrix 
X"^ = [xj^i\ij=i^,,,^n maps to a pair of u/rh paths {P2,Pi). Thus in the case that X = X'^ one 
has Pi = P2 so the mapping then is to a single family of paths. (The example of Figure ^ exhibits 
this.) Furthermore, for symmetric non-negative integer matrices X it is a known property of 
the RSK correspondence that pTj 

n n 

E^^-.:'=E(-1)''V„ (3.3) 
i=i i=i 



a fact which can also be derived within the setting of the PNG model (18|. Hence it follows that 



/i:^i</ 



[a"-2-gi...^J E «^^-^^"'^"'^^S/.(9i,...,9n) (3.4) 

/^:/^l<J 

where the second equality follows from ()2.7|) . Now we know from 7 that 

I n 

J2 a^U(-^y-'^^s,iqu . . . , g„) = ( n ((1 + + Qje''")) )^^^^. (3.5) 

IJ,:fil<l k=l j = l 

In 1)3. 5|) the average over the classical group 0(/) breaks into two parts, 

( • )o(0 = 2 " ~^ ^ ' ^*^"(')) 
where ( • )o+(z) denotes an average with respect to the eigenvalue p.d.f. for random matrices from 
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the classical group 0+(/), 



^;!J^!^nm->>m.j) n \^--e"'Cll-e'"'-">f, i even (3.6) 

^ ' ^' j=l l<j<k<l/2 



X JJ le". -e"'>|=|l-e'<'''+''»f , i odd, 

l<j<fc<(Z-l)/2 

(3.7) 

and ( • )o-{i) denotes an average with respect to the eigenvalue p.d.f. for random matrices from 
the classical group 0~(/), 

l/2-l 

X |e*^._e^e,|2|^_^.(^^,+e,)|2^ ^^^^^ 

l<j<'c<'/2-l 



(3. 



l<j<k<{l-l)/2 



(3.9) 



Substituting p.5() in p.4p we see that 



pv{ii^ <i)= ^K-''"]( ri(i ' (3-10) 

"■'"^ fc=l jr = l 



and substituting this in (|3.2|1 then gives 



I 

fe=i ^ ^ 

as first obtained by Baik and Rains [7j. 

As for the average (|2.5|) , the average (|3.11|) in the case q = also arises as a gap probability at 
the hard edge of a matrix ensemble. Thus consider the Laguerre symplectic ensemble, specified 
by the eigenvalue probability density function proportional to 

TV 

'=1 i<j<k<N 
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For a = 2{n — N) + 1, n > N, this is realized by the distinct eigenvalues of matrices of the form 
X^'X with X an n X N real quaternion Gaussian matrix (embedded as a complex matrix). Let 
E^{s;a;N) denote the probability that there are no eigenvalues in the interval (0, s). Then we 
know from |2U1 lllj that for a € Z>o 



^^'^<^(A;2a) := lim S[f4;2a;iV/2') =e-^/2/gV^E?iie''A (3.12) 

Af^oo VA^ / \ / 0{2a) 

thus coinciding with in the case a = and z'^ = X, I = 2a. 

3.3 The symmetry s 

Analogous to the case of the symmetry 0, let us suppose points are added below the diagonal 
from (0, 1) to (1,0) (to be referred to as the anti-diagonal) with rate zdz. Let the mirror images 
of these points in the anti-diagonal be also added, and let points on the anti-diagonal be added 
with rate f3dz. Then we know from (|3.1() and the following sentence that the probability there 
are exactly n points after time z, n — 2m of which are on the anti-diagonal, is equal to 

-Pz -z2/2^on-2m 

n\ 

Now the permutation of {1, 2, . . . , n} corresponding to a realization of this process is closely 
related to a realization of the symmetry 0. Thus let P = (^'(l), P{2), . . . , P{n)) be a permutation 
with the property that if P{i) = k then P{k) = j and thus correspond to the symmetry 0. 
Then P := {P{n),P{n - 1), . . . ,P{1)) has the property that if P{j) = k then P{n + 1 - k) = 
n + 1 — j and so corresponds to the symmetry s. Consequently the maximum length of the 
increasing subsequences in the case of s, say, is equal to the maximum length of the 
decreasing subsequences for 0. Furthermore, it follows from a theorem of Greene j24j (see 
also |35j ) relating row lengths of the standard tableau k corresponding to P, to increasing 
subsequences, and column lengths to decreasing subsequences, that the conjugate tableaux k' 
obtained by interchanging the row and column lengths in k, corresponds to P^. Also, it follows 
from (1231) that 



#(l's on the anti-diagonal) = y^(—l)-' ^kj (3.13) 
Setting k' = fi, these facts together imply 



fj,:fj,l<l 

P''~'"'qi---qn] P^^-'^-''>'"^'^s^{q,,...,qn) (3.14) 

fi:fil<.l 



(cf. 1)3. 4(1 ). In _7_ it was shown that 



21 



IJ.:fil<2l k=l ^ j = l ' 
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and 

n 21 n 

^t:^l<2«+l j = l k=lj = l 

where ( )sp(z) denotes an average with respect to the eigenvalue p.d.f. for random matrices from 
the classical group Sp(i), 



^ ' k=l l<j<k<l 



It follows immediately that 



Prills < 21) = J_[/3"-2-]/pj,^_^(^e-^ 

Sn.m. \f-^l- Pe ^l^* 



21 ^ ,21 



k=l j=l 



1 2i 



Pr(/^,„ < 2/ + 1) = + 



sp(0 



Thus with 



Pr(/^ < /) := e-'^-^Vs /3"-2'"s„,^Pr(/^,^ < /) (3.15) 



n=0 m=0 



we 



have [Z] 



21 

Pr(Z^<20 = e-/5-^V2/e.E?Lie-. TT 1 

\ 1 — ni 



_ /Je-^'^fe /spW 

Pr(;^<2/ + 1) = e-^'/2/g.E,tie'''A _ (3_lg) 

\ /sp{0 

There is an analogue of H3.12() for the second average in 1)3.16^ . For this consider the Laguerre 
orthogonal ensemble, specified by the eigenvalue probability density function proportional to 

AT 

l=\ i<j<k<N 

For a = {n — N — 1) /2, n > N, this is realized by the distinct eigenvalues of matrices of the form 
X^X with X an n X N real standard Gaussian matrix. Let -E'i'(s; a; N) denote the probability 
that there are no eigenvalues in the interval (0, s). Then we know from |2Ullllj that for a G Z>o 



A 

thus coinciding with (|3.16|) after setting X = z^, a = 1 



E\^"^{\-a) := Wni E\(—;a;N] =e-^/'^le^^f=^''^'') , (3.17) 

^ ^ ' N-*oo ^ \N J \ /Sp(a) ^ ^ 



3.4 The symmetry h 

For this symmetry only the points below both the diagonal and anti-diagonal are independent. 
Let us suppose points are added to this region with rate 2zdz. Suppose too that points are added 
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to the bottom half of the diagonal with rate adz, and to the bottom half of the anti-diagonal 
with rate pdz. The images of all these points must also be added. The probability of there 
being exactly 2n points after time z, 2m_|_ of which are on the diagonal, and 2m_ of which are 
on the anti-diagonal, is then equal to 



nl 

where 



7 9 

_f,-z -0!Z-f3z m+ am- j. 

I ce /J ''n,m_|_,m_ 



m-|-!m_!((n — m+ — m_)/2)! 
and it is required n — m+ — m_ be even. 

The permutation matrix P = [a;j.j]jj=i^...^2n corresponding to a realization of this process 
has the two symmetries P = P'^ = [xj^i]ij=i^,„^2n and P = P^ = [x2n+i-j,2n+i-i]i,j=i,...,2n- Now 
whereas a permutation matrix with the symmetry P = P^ maps under the RSK correspondence 
to a pair of identical standard tableaux (Ti , Ti ) say, a permutation matrix with the symmetry 
P = P^ maps to a pair of standard tableaux (Tg^, T2) where denotes the Schiitzenberger dual 
of T2 (see e.g. We note too that the number of permutation matrices of {1,2,..., 2n} with 

the symmetry P = P^ = P^ such that m+ (m_) members of {1,2, . . . ,n} have the property 
that P{i) = i {P{i) = 2n-|- 1 — i) is equal to tn,m+,m-- Hence with ln,m+,m- denoting the longest 
path length in a realization of the Hammersley process with symmetry ki, we have 

Pr(^.V,._ < I) = ^^f^ E a^?-(-)^-^-^7?^?-(-^)^-^'^^7^„ (3.18) 



where 



/2n ~ #(self dual standard tableaux, shape /i, content 2n). 
Analogous to (|2.6|) . let us define 

F/-((?i, . . . , g„) = E* qf' ■ ■ ■ "'^ (3.19) 

where the asterisk denotes that the sum is over all self-dual semi-standard tableaux of shape A 
and content 2n. (For self dual tableaux of content 2n, = #(2n + 1 — j)'s so there are n 
rather than 2n independent variables in the weightings.) Prom the definitions 

so we can write 



n,m_i_,m_ 



Baik and Rains jjj have provided the evaluation 

' 1 I id 

M:M<2Z j=l ^ k=l ^^^> 
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Hence 

If. iO ■ I 

Pvil^^. ^ < 21) = \a^^(r-\[ TT {^2zo^qX^. 

U{1) 



7=1 ^ k=l 



and consequently with 



oo 

Pr(/^<20:=e-^'— ^^^^ a"^+/3"^-tn,„^+,„._P<™^,^_ < 2^ 

we obtain the result |7j 

Fi{l^ < 21) = e-"'-"^-^^^ J] i:^^^e2^^fe=i=°'^^'=^^^^^. (3.20) 
Baik and Rains |7j have also provided the evaluation 

^l■.^lx<2l+l 

n I n 

k=l j=l k=l ' ' 

We readily deduce from this that [7j 

Pra^ < 2/ + 1) = Pr(/^ < 21) . (3.21) 

/3=0 

3.5 The symmetry b 

In the unit square suppose points are marked in the region y < 1/2 with Poisson rate 4zdz. 
For each point {x',y') so marked, also mark the image (1 — x',1 — y'), which corresponds to a 
reflection about the centre of the square. A realization of this process with 2n points corresponds 
to a permutation matrix with the symmetry X = (X^)^, or equivalently to a permutation of 
{1,2,..., 2n} with the property that if P{i) = j, then P{2n + l-i) = 2n + l- j. Note that 
there are (2n)!! = 2""n\ permutations of this type. 

Now we know [23 that if a permutation matrix P maps to a pair of standard tableaux 
(Ti,T2) under the Robinson-Schensted mapping, then P^ maps to (T2^,T/^) while P^ maps 
to (T2,Ti). Hence a permutation matrix with the symmetry X = (X^)'^ maps to a pair of 
standard tableaux {Ti,T2) constrained so that 

Ti = T^, T2 = T2^. 

Consequently with denoting the longest path in a realization of 2n points, 

fi.fj,i <l 

= -^^[ai ■ ■ ■ anbi ■ ■ ■ bn] E ^^("i; • • • > a'„)s^(6i, . . . , 6n)- 

fiifii <i 



n=0 r?i+,m_>0 
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Baik and Rains [7j have derived the results 



L5 • • • ) bn) 



s^(ai, . . . ,a„)s^(6i, 

fi:fil<2l 

n I 

1111(1 + a,e^'-)il + b,e~ 
j=i k=i 



n 21 



U{1) 



j=i k=i 



U{l)®U{l) 



s^{ai,...,an)s^{bi,...,hn) 

fi:fMl<2l+l 

n l+l 

lll[{l + aje''-){l + b,e-''''] 
j=i k=i 

n 21+1 



n I 



U{l+l) 



\{\[{l + a,e''-){l + h^e-^'-) 



j=i k=i 



1111(1 + a,e^'^){l + b,e-^'-: 



j=i k=i 



c/(/+i)ec/(/) 



Consequently 



and thus with 



we see that 



Piil^^ < 2/ + 1) 



-(\y 

"n! \l ^ 

k=l 
, 21+1 

—(\y 



2n 



2n 



k=l 



,2 ^ 2"^2n 



U(l)(BU{l) 



U{l+l)®U{l) 



Pr(/- < /) := ^ ^P^(^2n < 



n=0 



Pr(/°<20 = [Vv{P <l) 



Pr(F < 2/) = Pr(F < / + 1) Pr(/° < /) 



where Pr(P < I) is given by H2.5|) . 



U{1) 



(3.22) 



(3.23) 



4 The Hammersley process with sources on the boundary 

The Hammersley process which relates to the scaling function for one-dimensional stationary 
KPZ growth is the original model (Poisson points in a square), generalized to allow independent 
Poisson rates for points forming on the boundaries y = (with intensity a_) and x = (with 
intensity a+). In the PNG model picture, these boundary points correspond to growth from 
the boundary of the expanding droplet, with rate a^dt on the left boundary, and rate a-dt 
on the right boundary. Such boundary growths can be realized by an initial condition of mean 
slope — a_ («+) for X < (x > 0), created by a staircase structure (downward sloping for 
X < 0, upward sloping for x > 0), with vertical increments of one unit with intensity a_ (a+). 
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A formula for the cumulative probability of the longest increasing subsequence length of this 
model has been obtained by Baik and Rains This was obtained as a limiting case of an 
inhomogeneous lattice generalization of the Hammersley process due to Johansson j29j . 

In the Johansson model, on each site < i,j < n of a (n + 1) x (n + 1) square grid 

there is a non-negative integer geometric random variable Xij with parameter aibj so that 

Pr(x,j =k) = {l- a,hj){a,hjf. (4.1) 

A type of directed last passage percolation is to form an u/rh lattice path from the site (0,0) 
to the site (n, n) such that it maximizes the sum total of the random variables associated with 
the sites. Of interest then is the value of this maximized sum, 

L(n,n):= max Xjj. 

(0,0)u/rh{n,n) 

It is not hard to show that in the RSK correspondence mapping the non-negative integer matrix 
X = [xij]jj=o,...,n to a pair of semi-standard tableaux of the same shape and content n + 1, the 
length of the first row of the tableaux is equal to L(n, n) (one approach is to show that L{n,n) 
and the height at the origin in the PNG model picture of the RSK correspondence satisfy the 
same recurrence). Further, one has that with entries of the matrix chosen according to (|4.1() . 
the RSK correspondence maps to pairs of weighted semi-standard tableaux, weights {aj}, {bj} 
repsectively. Consequently 

n 

Pr(L(n,n) </) = J] (1 - a,6fc) 

j,k=0 t^:tJ.i<l 

and thus according to (|2.8jl 

n n I 

Pr(L(n, n) < = n (1 - «^-^^)( 11 11(1 + a,e'^'=)(l + b,e-^'')) . (4.2) 

j,k=0 j=0 k=l ^ ' 

Suppose we set 

ai = bi = t/n (i = l,...,n) (4.3) 

and take the limit n ^ oo. Then the probability that there are k points in the region 1 < i,j < n 
has the large n behaviour 

where the first factor corresponds to the number of different ways of choosing k sites from the 
grid of sites, the second factor is the probability that those k sites are occupied, while the 
final factor is the probability that the remaining sites are empty. Hence we reclaim the setting 
of the original Hammersley model with Poisson parameter t^. Indeed with ao = 6o = we see 
that by substituting H4.3() in (|4.2() and taking the limit n ^ oo H2.5|) results with X = t. 
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Instead of setting ao and bo to zero, suppose we choose oq = a+, = ai- Then the argument 
leading to (|4.4j) shows that along x > (y > 0) we obtain a Poisson process of intensity a+t 
(a^t). At the origin there remains a non-negative integer variable chosen according to the 
geometric distribution with parameter a+a_, g{a+a-). Let us denote by L"^(t, a+,a_) the 
length of the longest up/right path in this process, and let us denote by a+,a_) the same 
quantity but with the geometric random variable at the origin removed. Clearly 

L+(t,a+,a_) = L(t,a+,a_) +x, X G 5(a+a-), (4.5) 

while the appropriate limit of 1)4.2(1 gives 

Pr(L+(t, a+, a_) < /) = (1 - a+a_)e-(°++"-)*-*' A (4.6) 

where 

I 

Di = (l[{l + a+e'''^){l + a-e~'^^)e^'^U'°'^A . (4.7) 
To now obtain a formula for Pr(L(t) < I), introduce the generating functions 

oo oo 

Q(x) = Pr(L(t) < l)x\ g+(x) = Pr(^+(i) < l)x^- 
1=0 1=0 

Then using (|4.5() we see that 

oo I 

Q+(x) = ^x'^Pr(L(t) < / - A;)Pr(x = A;) 
1=0 k=0 

oo I ^ 

= (1 - a+a_) y y Pr(L(t) < I - k){a^a.f = ~ Q{x) 

1=0 k=0 ^ 

where the final equality follows by writing x' = x^^^x^ and summing independently over I — k 
and k. Multiplying both sides of this identity by 1 — xa+O- and equating like powers of x gives 

M 

Pr(L(t, a+,a_) < /) = e-(°++"-)*-*' (A - a+a_A_i). (4.8) 
To proceed further, Di is expressed in terms of 



Di = Di 



a+=a-=0 \ / U{1) 



and monic orthogonal polynomials {tTj (e*^)}j=o,i,... with respect to the weight e^*™''^, 

]- r 7r,(e^^)^^^e2*™«^ = ^5,- fc. (4.10) 

''''' J -TT 



2ir 

For this purpose, let 



Kiz) := z"7r„(z-i). (4.11) 
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Then with {7rj(e*^)}j=o,i,... the monic orthogonal polynomials corresponding to the general real 
weight w{6) replacing g^*™*^^ in (|4.1()j) . we have the general formula 



^^■'-^ ^ ^'^ ''u(n) _ <+i(a;X+i(y) -7r„+i(x)7r„+i(y) 



^ <(^X(j/) - Xy'Knix)'Kn{y) 

1 — xy 

1 " 

— Y.'^l'^k{x)7rk{y). (4.12) 



" fc=0 



Application of the second of these equalities shows 



= M 

1 — a^a_ 



and in particular, applying I'Hopitals rule, we see that 



= \il-l)nii-a)nii-a~^)-an'ii-a)TTii-a-^)-a-^nii-a)7rl{-a~^)}Di. (4.13) 
i_i_=i/«_=« J 



The scaling behaviour of Di = Di{t) as defined by (|4.9I) is the result ()2.1U() of Baik, Deift 
and Johansson. Thus we have 

e"*'z)[2,+,V3,](t) =i^GUE(s). (4.14) 
We note that -Fgue('S) has the exact evaluation [HTj 

i^GUE(s) =exp{ {t-s)q^{t)dt^ (4.15) 

where q{t) is the solution of the non-linear equation 

= tq + 2q^ (4.16) 

(a special case of the Painleve II equation) subject to the boundary condition 

q{t) ~ -Ai(t) as t ^ oo (4.17) 

where Ai(t) denotes the Airy function. 
In addition to setting 

l = [2t + t^/^s], (4.18) 
to obtain critical scaling behaviour, the parameter a must be related to i by [S] 

a = 1 - y/t^/^ (4.19) 

where y is fixed. To specify the corresponding scaled form of the orthogonal polynomials in 
(j4.13|) . one establishes equations for their variation with respect to the scaled variables s and y 
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jnilSES]- In relation to the former, one recalls that in general monic orthogonal polynomials on 
the unit circle satisfy the coupled recurrences (SO] 

iTn+i{z) = Z7r„(z) + r„+i7r*(2;) 

<+l(^) = rn+lZTTn{z) + Tt'^{z) (4.20) 

where r„ = vr„,(0). For the weight w{z) = g^*™*^^ the scaled form of r; = ri{t) with the 
substitution ()4.18|) can be determined from a difference equation — a form of the discrete 
Painleve II equation — satisfied by {r„} jSH [23 E3 |2H1 H US IH EU • 

Proposition 1. With r.„ := 7r„(0), the sequence {r„} for the polynomials with orthogonality 
i4-l(^ satisfies a form of the discrete Painleve II equation 

n r„ 



tl-ri 

subject to the initial conditions 



rn+i + rn-1 (4.21) 



^0 = 1, n = -y^, (4.22) 

where Ii,{t) denotes the Bess el function of pure imaginary argument. 

Proof. Following |2H| j let us show how 1)4. 21() can be derived using simple properties of the 
weight 

w{z) = e'^'+^/'\ (4.23) 

together with a general formula from the theory of orthogonal polynomials on the unit cir- 
cle. Regarding the latter, set (j)j{z) := KjTTj{z) so that according to ()4.10|) {(j)j{z)} defines an 
orthonormal set of polynomials on the unit circle, and introduce the coefficient by 

= + + • • • + cPniO). (4.24) 

Then it is generally true that [36j 



^=gr..i^„ rr.= ^. (4.25) 



j=0 J 

Specific to the weight (|4.23|1 consider 

d , ,\ , , ,- ^ dz 



J 



z'^(^ — w{z)j(l)n{z)(t)n+l{z) 



2TTiz 

where C is a simple closed contour encircling the origin. Noting from (|4.23|) that 

^w{z)=t[l-^)w{z) 
we see from integration by parts, the structure 1)4. 24() . and the orthonormality of {(j)j{z)} that 

j = _(n + l)^ + (n + l)^. (4.26) 
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On the other hand, direct evaluation of J using H4.24() and the orthonormahty of {(j)j{z)} shows 

j = t(_f^_i!^J!!^). (4.27) 



Equating (|4.26|) and (|4.27|) and ehminating /„ using (|4.25jl (with rj = rj) we arrive at (|4.21j) . □ 
Following ^37, with 

Rn{t) := 

and making the ansatz 

%+tV3,](t) ~ r^/Vs), t^oo, (4.28) 



we see that formally the difference equation (|4.21j) becomes the differential equation 

--r = su + 2u^. (4.29) 
Further, for 1)4. 28(1 to be compatible with the first of the initial conditions ()4.22|1 . one must have 



u{s) ~ -J-s/2. (4.30) 

Now (|4.29() is the same particular Painleve II equation as (|4.16|) . In fact it is a celebrated 
result |26j in the theory of the Painleve II equation that 1)4. 16() has a unique solution with the 
asymptotic s — oo behaviour (|4.3()|1 . and the asymptotic s — > oo behaviour u{s) ~ — Ai(s) as 



in (|4.17j) . Thus we conclude 

^[2t+tV3,](t) ~ t-'/'g(s), t^oo. (4.31) 

According to (|4.13j) . our interest is in 7r„(— a) and 7r*(— a), so we should set z = a in (|4.2Uj) . 
Introducing 

Pn{a) = e-*°<(-a), Q„(a) = -e-*"(-l)"7r„(-a), (4.32) 
one sees that (|4.2U() is consistent with the existence of the scaled quantities 

a{s,y) := ^Inn (1 " v/t'^'), b{s,y) := ^hm Q[2,+,V3,] (1 - v/t'^'), (4.33) 
and that furthermore (|4.2U|) reduces to the partial differential equations 

da , db , / , . X 

IT = TT = 9" - 4.34 
OS us 

where use has also been made of 1)4. 28p and 1)4. 31(1 . We note too that existence of the limits 
(|4.33|1 together with the formula 

1 J 



/Di = ^ y^K^7rfc(a+)7rfe(l/a_ 



a+=l/a.=a — 

which follows from the final equality in (|4.12j) . allows for the formal derivation of the limit 
g(s,y):= lime-(-++--)*5[2i+iV3,](t) ^ ,^,JD^2t+t^/s^^ 

a{s',y)a{s',-y)ds' = b{s' ,y)b{s' , -y) ds' (4.35) 
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Similarly, at the same formal level, we see from H4.8() . the definition of g{s,y), and ()4.14|1 that 



Fy{s) := ^hmPr( -^^ <. 



= —(^g{s,y)FGVEis)). (4.36) 

To fully determine a and b, and thus g{s,y) and the scaled distribution Fy{s), it remains to 
specify equations for their dependence on the scaled variable y as introduced in (|4.19|) . Such 
equations follow from differential equations in z for 7r„ and vr* [281 IH EH] . 

Proposition 2. We have 

7r„(z) = [- + ^ )7^n(z) + [— ^j^n(^) (4-37) 

K'iz) = (-^+r„t)vr„(.) + (-t+^i^)<(z). (4.38) 

Proof. Let us show how ()4.37() can be derived using results from j28| . In terms of the 
polynomials (|mH) for the weight it is shown in [251 eq. (2.79) with t 2t] 



, / , ^ Kn-1/ , t K„_i</>„_i(0) >„+i(0)</>„(0)\ 
(AnA^) = [n+- + t — — t )4>n-l{z) 

t Kn-1 0„_l(O) 

0n(z). (4.39) 



Z Kn (/>n(0) 

But for a general weight |36j 

K;n-l2</'n-l(^;) = Kn<pn{z) - 0„(O)0* (z), 1 - = ('^lizi') , 

which together with show (|0^ reduces to (|07)) . □ 

The differential equations 1)4. 37|) , 1)4. 38() are consistent with the existence of the limits ()4.33() , 
and furthermore assume the scaled form 

—a = q^a - {q' + yq)b, —b = {q' - yq)a + (y^ - s - q^)b. (4.40) 

Together (|4.34|) and (|4.4()|) determine the scaled quantities (|4.33j) once appropriate initial con- 
ditions are specified. 

According to 1)4. 19() . when y = 0, a = 1, so according to ()4.33|1 we require the behaviour of 
7r„(-l), <(-!). First we note from (HTH) that <(-!) = (-l)''7r„,(-l) and thus (lOl gives 
a(s,0) = -6(s,0). It then follows from that 



a{s,Q) = Ae-^^'\ U{s) = -J q{t)dt. (4.41) 

To determine A, we note that formulas in imply 

lim (-l)«vr„(-l) = e* 
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so we see from 1)4. 32() and (|4.33() that 



a{s,0) — > 1 as s — > oo 

and so 

a{s,0) = -b{s,0) = e~^^'\ (4.42) 

The quantities o(s,y), b{s,y) are now fully determined, and thus so to is g{s,y) as specified by 
()4.35|) and the scaled cumulative distribution Fy{s) as specified by ()4.36|) . 

The scaled quantities satisfy a number of further properties of interest. First, with the initial 
condition 1)4. 42() it is easy to see from ()4.34|1 and ()4.4U|1 that a(s, y) and b{s, y) are related by 

a{s,y) = -b{s,-y)ey"-'y. (4.43) 

From this it is simple to verify that ()4.35|) can alternatively be written [7] 

d d 
9{s, y) = a{s, -y)—a{s, y) - b{s, -y)—b{s, y), (4.44) 

which relates to Second, it follows from (jUTZj), and that 

b{s,y) ~ -e~y^. 

s—^oa 

One then sees from (|4.43|) . (|4.4fl|) and (|4.44|) that 

g{s,y) ~ (4.45) 

Using this, the mean of the distribution dFy/ds can be computed as 

j sdFy = _^liin {sFy{s) - g{s, y)FGUE(s)) = 



where the first equality follows from integration by parts and 1)4. 36() . while the second equality 
follows from the facts that Fy[s) and i*GUE(s) approach 1 exponentially fast, together with 
(|4.45|) . This then motivates defining the shifted cumulative distibution |33j 

Fy{s) ■.= Fy{s + y^) 

for which the corresponding distribution function dFy/ds has mean zero. From the discussion 
in [HH], ones sees that of immediate interest to KPZ growth is the second moment 

/oo 
S^dFy, 
-00 

which is a function of the scaled parameter y. Prahofer and Spohn have used the results 
revised above, evaluated using high precision computing, to accurately tabulate this quantity 
and discuss its properties. 
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5 Scaled limits for the symmetrized models 



In a remarkable analysis Baik and Rains |S] have provided similar evaluations to the one detailed 
in Section 4 of the scaled limits of the cumulative distributions in Section 3. Here we will be 
content with drawing attention to a subcase of two of these: Pr(Z^ < I) in the case a = and 
Pr(/^ < 21) in the case P = 0. Thus it was first proved in [S] (see for a subsequent simplified 
derivation based on the identities (I3.12I) and (|3.17j) ) that 



Here Fqse{s) denotes the scaled cumulative distribution of the largest eigenvalue for large Her- 
mitian random matrices with real quaternion elements, while i*GOE(s) denotes the same for large 
real symmetric matrices. 

Both -FgseI-s) and Fqoe{s) can be expressed in terms of the same Painleve 11 transcendent 
as in the evaluation H4.15() of -FgueI-s)- Thus one has jSHI 



The origin of such an identity can be traced back to a special property of a particular marginal 
distribution of the joint probability for the row lengths of semi-standard tableaux relating to the 
symmetry s j3 1161 118j . The marginal distribution is defined by summing over every second row 
of the semi-standard tableaux. We also draw attention to the fact that Fgue(s) and -Fgoe(s) are 
T-functions for certain Painleve II systems j37l I19j . Similarly, -Fgse(s) is the arithmetic mean 
of two r-functions, both of which correspond to Hamiltonians satisfying the same differential 
equation, differing only in the boundary condition [19| . 
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